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ABSTRACT 



We present SuperWASP observations of HAT-P-14b, a hot Jupiter discovered by Torres et 
al. The planet was found independently by the SuperWASP team and named WASP-27b after 
follow-up observations had secured the discovery, but prior to the publication by Torres et al. 
Our analysis of HAT-P-14/WASP-27 is in good agreement with the values found by Torres et al. 
and we refine the parameters by combining our datasets. We also provide additional evidence 
against astronomical false positives. Due to the brightness of the host star, V mSLg = 10, HAT- 
P-14 is an attractive candidate for further characterisation observations. The planet has a high 
impact parameter, b = 0.907 ± 0.004, and the primary transit is close to grazing. This could 
readily reveal small deviations in the orbital parameters indicating the presence of a third body 
in the system, which may be causing the small but significant orbital eccentricity, e = 0.095 
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± 0.011. The system geometry suggests that the planet narrowly fails to undergo a secondary 
eclipse. However, even a non-detection would tightly constrain the system parameters. 

Subject headings: planetary systems — stars: individual: (HAT-P-14, WASP-27, GSC 3086- 
00152) - techniques: spectroscopic, photometric 



1. Introduction 



There has been a rapid increase in the number of transiting pla nets disc overed each year d ue to 
dedicated ground- and space- b ased surveys: HAT (IBakos et al.ll2002r). TrES (Alonso et al. 12004). XO 
JlVlcCullough et al.l2005h . WASP JPollacco et al.l2006h . CoRoT jBaglin et al.l2006h and Kepler jBorucki et al 
201ol) . This trend looks set to continue, with the discovery of over 35 new planets published already this year 



(mid 2010), which represents more than a third of the total number of transiting planets known. We have 
now surpassed one hundred transiting planet discoveries. By increasing the number of well-characterised 
transiting planets, we are beginning to be able to test theories of planet formation and evolution on a secure 
statistical basis. For example, the obliquity of hot Jupiters appears t o show that the mechanism by which 
plane ts migrate cannot be explained solely by planet-disc interactions ( Fabry cky & Winnll2009 : Triaud et al. 



2010). 



Transiting planets orbiting the brightest stars are extremely valuable due to the extensive follow-up ob- 
servations which can be undertaken. This includes secondary eclipse observations to determine atmospheric 
properties and refine orbital parameters. In order to find planets orbiting bright stars, a wide field-of-view 
must be employed to encompass enough stars to make the detection probability significant. It is therefore 
inevitable that some areas of the sky will be monitored by multiple groups, and indeed se veral transit- 
ing p lanets have been discovered independently: WASP- lib ([West et al.l l2009) = HAT-P-lOb (|Bakos et al 



2009) and XO-5 (IBurke et al.ll2008l : IPal et al. 2009). 



Here we describe observations of the transiting planet HAT-P-14b (ITorres et al.ll2010bl . hereafter T10) 
by the SuperWASP survey. The star was monitored by the SuperWASP station on La Palma and followed up 
with photometric and spectroscopic observations using the Faulkes Telescope North, Liverpool Telescope, 
Observatoire de Haute-Province and Nordic Optical Telescope. The planet was named WASP-27b after the 
detection was secured, prior to the publication of the T10 paper. While the HAT-P-14b designation is to be 
preferred in recognition of priority of publication, it is important that we retain WASP-27b as a secondary 
designation for statistical completeness in studies of the physical properties of independently-discovered 
WASP planets. 

Section |2] describes the observations of HAT-P-14b/WASP27b, including the SuperWASP discovery 
data and photometric and spectroscopic follow-up. The results of the derived system parameters are pre- 
sented in Section [3l including an independent spectral analysis and a joint fit with the T10 data in order to 
refine the planetary parameters. We discuss our findings and compare them to those found by T10 in Section 

a 
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2. Observations 



2.1. SuperWASP photometry 

HAT-P-14 is a bright, (V mag = 9.98) star in the constellation Virgo, located at a J2 ooo = 17 h 20 m 27 s .87, 
feooo = +38°14 / 31".9 (GSC 3086-00152; 2MASS 17202788+3814317). It has been observed by Super- 
WASP (La Palma) since 2004 and a total of 25,474 photometric points have been collected. A periodic 
signature was detected in the observations obtained betwe en March and August in 200 7, 2008 and 2009. All 
the data were processed using the pipeline described in ICollier Cameron et al.1 (120061) . yielding a detection 
of a transit-like feature with a period of 4.6278 d in multiple cameras and successive seasons. The folded 
light curve is shown in Figure Q] 



The target underwent several consistency tests aimed at eliminating false positives (see lCameron et al. 



20071) . Resolved blends were ruled out by inspecting Digitized Sky Survey (DSS) images which showed 
the star to be isolated within the WASP photometric aperture. No significant ellipsoidal variability was 
measurable in the WASP light curve folded on the transit period. The transit depth and duration yielded a 
planet-like radius for the companion, and a stellar density appropriate to a main-sequence host star of the 
effective temperature derived from the 2MASS colours. Having passed all these tests, the star was selected 
for follow-up observations. 



2.2. Photometric follow-up 

We obtained six further light curves in order to refine the photometric parameters, shown in Figure [2] 
All photometric data presented in this paper are available from the NStED database^ 



http://nsted.ipac.caltech.edu 



o 
E 




0.3 



Fig. 1. — The combined, unbinned SuperWASP light curve for HAT-P-14, folded on the orbital period of P 
= 4.627654 d. Superimposed is the model transit light curve, based on the system parameters determined 
from a global fit, see Section [3721 
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Fig. 2. — Photometry of transit events of HAT-P-14b taken with FTN, Astrograph and LT in various bands. 
The light curves have been offset from zero by arbitrary amounts for clarity. Superimposed are the model 
transit light curves based on the determined system parameters (see Section l3^2l . 
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Two full transits were observed with Faulkes-North (FTN, Hawaii), on 2009 September 18 (SDSS 
i-band) and 2010 May 21 (Pan-STARRS-z filter). The first transit was taken using the Merope camera 
(EM01) in binning 2x2 mode and defocussed to enable an exposure time of 40s. Aperture photometry 
was performed using an 18-pixel radius aperture and ensemble photometry using the 2 bright, non- variable 
comparison stars within the field-of-view. By 2010 May 21, the new Spectral camera (FS03) was available at 
FTN and defocussed to obtain 20s exposures. The 10.36 x 10.22 arc-minute field-of-view included 4 bright 
comparison stars. These observations took advantage of the 4ag autoguider to maintain the locations of the 
stars to within an RMS ~ 2 pixel of the same location throughout the observations and led to a noticable 
improvement in the photometric precision. All the FTN data were reduced using the pipeline written at 
Liverpool John Moores University, and photometry performed using the IRAF/DAOPHOT package. 

A partial transit on 2010 March 13 was observed simultaneously by the robotic 2.0-m Liverpool Tele- 
scope (LT) and an 1 8cm Takahashi Astrograph telescope, both located on La Palma. Two further full transits 
were observed, on 2010 May 2 using the Astrograph, and on 2010 June 22 by the LT. 

We used an Andor 1 x IK e2v CCD detector with a 5 s readout time mounted on the Astrograph. For 
the first transit, the telescope was defocussed with FWHM = 20" and 685 exposures of duration 20 s were 
obtained in the 4.8 hour period. For the second transit, the exposure time was increased to 30 s and FWHM 
= 28"giving 456 observations in the 4.5 hour period. The data were reduced and differential photometry 
performed using the PYRAF/DAOPHOT package. The observations were subsequently binned with 3 and 
2 points per bin, respectively. 



We also used the RISE frame transfer CCD located on the LT dGibson et alj|2008l : iSteele et al.1 12008). 
For both transit observations, the telescope was defocused by -1.2mm giving FWHM =18 and 21 pixels on 
the respective nights (0.54 arcsec/pixel). The CCD was binned 2x2 and exposures of 3.7 s taken with effec- 
tively no dead time. During the 3.0 h and 2.3 h periods, 3520 and 3290 frames were taken, res pectively. The 



data were reduced and differential photometry performed using the ULTRACAM pipeline (|Dhillon et al 
12007m . and then binned with 10 points per bin. 



2.3. Spectroscopic Observations 

We obtained follow-up spectroscopic observations to determine the planetary, orbital and stellar pa- 
rameters. The star was initially observed using the Fibre-fed Echelle Spectrograph (FIES) mounted on the 
2.5-m Nordic Optical Telescope. In total, 24 usable spectra were obtained between 2009 July 12 and 2010 
April 05. Two spectra were removed from the analysis due to Moon contamination affecting the radial 
velocity determination. FIES was used in medium resolution mode (R = 46,000) with simultaneous ThAr 
calibration, and observations were conducted using exposure times of 1200s. We used the bespoke data 
reduction package FIEStoo|§to extract the spectra. An IDL cross- correlation routine was used to obtain ra- 
dial velocities (RVs) by fitting gaussians to the cross-correlation functions (CCFs) of 30 spectral orders and 



2 http://www.not.iac.es/instruments/fies/fiestool/FIEStool.html 
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Table 1: Radial velocities (RV) and line bisector spans (V span ) of HAT-P-14. 



BJD 


RV 


Error 


T r 

v span 


Instrument 


-2 450 000.0 


(kms ) 


(kms ) 


(kms ) 




5024.5363 


-20.6620 


0.0071 


-0.0031 


FIES 


5025.5236 


-20.4951 


0.0104 


0.0030 


FIES 


5026.4324 


-20.2582 


0.0105 


-0.0155 


FIES 


5026.6481 


-20.2427 


0.0145 


0.0010 


FIES 


5040.5041 


-20.1997 


0.0179 


-0.0326 


FIES 


5041.5629 


-20.3024 


0.0178 


-0.0260 


FIES 


5085.3553 


-20.5455 


0.0072 


0.0119 


FIES 


5086.3534 


-20.2943 


0.0150 


-0.0064 


FIES 


rAfin <~i c f\c\ 

5087.3509 


-20.2642 


0.0126 


-0.0072 


FIES 


5096.4062 


-20.2160 


0.0147 


0.0046 


FIES 


5098.3406 


-20.6554 


0.0073 


-0.0336 


FIES 


C f\C\C\ O O /" 1 

5099.3361 


-20.5338 


0.0088 


-0.0087 


FIES 


5100.3558 


-20.2565 


0.0127 


-0.0021 


FIES 


5101.3677 


-20.2778 


0.0121 


-0.0029 


FIES 


5H9.3355 


-20.2386 


0.0150 


-0.0014 


FIES 


5291.6152 


-20.3542 


0.0094 


-0.0026 


FIES 


5291.6300 


-20.3725 


0.0080 


0.0090 


FIES 


/'A AC\ 

5291.6449 


-20.3932 


0.0075 


0.0016 


FIES 


5291.6597 


-20.3904 


0.0095 


r\ r\r\i—i a 

-0.0074 


FIES 


5291.6146 


-20.3904 


0.0063 


-0.0294 


FIES 


5291.6895 


-20.3837 


0.0079 


-0.0231 


FIES 


5291.7043 


-20.4016 


0.0086 


-0.0029 


FIES 


jzy i. / iyz 


-ZU. DO J I 


A AA^T 

U.UUuJ 


A AAC/I 




5291.7340 


-20.3984 


0.0065 


0.0062 


FIES 


5068.3272 


-20.1678 


0.0146 


0.0236 


SOPHIE 


5071.4331 


-20.5014 


0.0123 


0.0540 


SOPHIE 


5072.3251 


-20.2567 


0.0126 


-0.0163 


SOPHIE 


5073.3273 


-20.1498 


0.0166 


-0.0242 


SOPHIE 


5074.4166 


-20.4524 


0.0133 


0.0372 


SOPHIE 


5075.4225 


-20.6275 


0.0132 


0.0031 


SOPHIE 


5076.3691 


-20.4281 


0.0192 


-0.0165 


SOPHIE 


5079.3532 


-20.5479 


0.0144 


0.0437 


SOPHIE 


5080.3964 


-20.5847 


0.0162 


0.0308 


SOPHIE 
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Fig. 3. — Upper panel: Phase folded radial velocity measurements of HAT-P-14, combining data from 
NOT/FIES (open squares) and OHP/SOPHIE (filled circles), and superimposed with the best-fit model RV 
curve based on the determined system parameters. The systematic velocity for each dataset was subtracted 
from the RVs. Middle panel: Residuals from the orbital fit plotted against time. No long-term trend is visi- 
ble. Lower panel: The bisector span measurements as a function of radial velocity, showing no correlation. 
The uncertainties in the bisector spans was taken as twice the RV uncertainties. 
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taking the mean. A template spectrum was constructed by shifting and co-adding the spectra, against which 
the individual spectra were cross-correlated to obtain the final velocities. The template was cross-correlated 
with a high signal-to-noise spectrum of the Sun to obtain the absolute velocity to which the relative RVs 
were shifted. The RV uncertainty is given by a = RMS(u) / vN, where v is the RV of the individual or- 
ders and N is the number of orders. Line bisector spans (V Bpan ) were computed using the difference in the 
position of the midpoint of the CCF at 25% and 75% of its depth. 

A further nine spectra were taken between 2009 August 24 and 2009 September 05 wi th the stabilised 



echel l e spectrograph SOP HIE at the 1.93-m telescope of Observatiore de Haute-Provence (IPerruchot et al. 



2008: iBouchy et all 120090 . The spectrograph was used in high efficiency mode (resolution R = 40,000) 



and the observations w ere taken with a sig nal-to-noise ratio S/N~30 to minimise the Charge Transfer In- 
efficiency (CTI) effect (|Bouchy et al.ll2009l) . Two 3 arc-second diameter optical fibres were used, the first 
centred on the target and the second on the sky to simultaneously measure the background to r emove contam 



ination from scattered moonlight. The spectra were extracted using the SOPHIE pipeline (IPerruchot et al. 



2008) and RVs computed from a weighted cro s s corr elati on of each spectru m with a numerical mask of 



G2 spectral type, as described by iBaranne et all (| 19961) and lPepe et all (|2002h . Cross correlation with other 
masks (F0, K5 and M5) produced an identical RV semi-amplitude, indicating that the variation is unlikely 
to be caused by a blended eclipsing binary system of unequal masses. 

The radial velocities and line bisector spans are shown in Table Q] and plotted in Figure [3] A possible 
correlation between the bisector span and radial velocity was noted by T10 but was attributed to moonlight 
contamination rather than to a blended eclipsing binary. We find no correlation, supporting the signal's 
origin as a planetary companion rather than an astronomical false positive. 



2.4. High resolution imaging 

In order to search for additional close stars which may cause blends, we obtained high-resolution H- 
and A'-band imaging with the near-infrared camera INGRID, fed by the adaptive optics system NAOMI on 
the 4.2-m William Herschel Telescope (WHT), see Figure HI Images were taken in two sky position angles 
(PA, 0° and 30°) in natural seeing of 0.5", and the corrected PSFs have a FWHM of 0.12"and circular 
contours. No contaminating sources are detected to limiting magnitudes of H ~ 17.1 ± 0.25 and K s ~ 16.7 
± 0.25 at distances > 4 x FWHM from the centre of the corrected profile. At 1.5 x FWHM the limiting 
magnitude is 2 magnitudes brighter than this. 



3. Results 

3.1. Stellar parameters 

We derived the fundamental stellar parameters of HAT-P-14 using the following techniques, and the 
results are shown in Table |2] First, we performed a spectral analysis using SME (Spectroscopy Made Easy, 
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Table 2: Stellar parameters for HAT-P-14. 



Parameter (Unit) 



Value 



Photometric & spatial properties: 

RA (J2000) 17:20:27.88 

DEC(J2000) +38:14:31.7 

V(mag) 9.98 ± 0.058 

J (mag) 9.094 ± 0.021 

H (mag) 8.927 ± 0.020 

K s {ma.g) 8.851 ±0.019 

MRA (mas year -1 ) 1.6 ±0.7 

/"dec ( mas year -1 ) -4.7 ± 0.7 

[/(kms -1 ) 6.5 ±1.0 

^(kms- 1 ) -10.6 ±0.1 

^(kms- 1 ) -6.0 ±0.8 

Galactic longitude (deg) 62.6 

Galactic latitude (deg) 33.5 



Spectroscopic properties: 

T cS (K) 
log 5 (cgs) 
v mic (kms -1 ) 

"mac (kms -1 ) 

vsini (kms -1 ) 
[M/H] 



Derived properties: 

M* (M ) a 

R* (R ) b 

P* (Pq) c 
U (L Q ) d 
Age (Gyr) d 
Spectra type 



6583 ± 100 
4.1 ±0.12 

0.85 

5.28 
8.4 ± 1.0 
0.08 ±0.1 



1.30 ± 0.03 
1.48 ± 0.05 
0.398lffi 



3.91 
1.6 
F5V 



0.37 
+0.54 
+0.4 
0.3 



Note: The photometric and sp atial properties are t aken from or derived using the following sources: V 
(TASS.IPJchmond et al.l l200oh. J, H, K s QMASS. ISkrutskie et alJbood). proper motions (NO MAD. 



Zacharias et al. 



2005J), space velocities (calculated according to ljohnson & Soderblomll 19871. updated to 



J2000, for a right-hand c oordinate system and corrected for solar motion using data from 



Dehnen & Binneylll998l) and galactic co-ord inates (NED ). 



a derived from the empirical relationship of iTorres et al.1 (12010a!) . 

b determined using M* and p+. 

c found using the lightcurve geometry of WASP data. 

d from stellar models 
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see IValenti & Piskunovlll996r). followi ng the method of IValenti & Fischer! (120051) and using the Kurucz 
model atmospheres ( Kurucz et alll984 ). Individual NOT spectra were normalised and co-added to produce 
a single high signal-to-noise spectrum. The H a , Na I D and Mg I b lines were fitted simultaneously to 
determine the spectral parameters. This yielded the values T e ff = 6583 ± 100 K and logg spcc = 4.02 ± 0.1, 
which indicate that the star is of spectral type F5V (IGrayl 120081) . Values for microturbul ence (v m \ r . = 0-85 
kms" 1 ) and macroturbulence (t> ma c = 5.28 kms -1 ) are taken from the prescription of IValenti & Fischer 
(|2005r) . The projected stellar rotation velocity (v sin i) was determined to be v sin i = 8.4 ± 1.0 km s _1 and 
the metallically [M/H] = 0.08 ± 0.10. 

The mean stellar density (p*), found from the light curve geometry, can be used as a luminosity indi- 
cator for stellar e volutionary models and often provides a stronger constraint than the value of log g from 
spectral analysis (|Sozzetti et alJl2007T) . We used a Markov Chain Monte Carlo (MCMC) approach to glob- 
ally model the photometric and radial velocity data (see Section I3T21 and obtained a mean stellar density of 
OJQS^nq ^ 37 f0- W e then compared T c g, [M/H] and p+ with the theoretical stellar evolutionary models of 

1.6i^Gyr, 



Girardi et al.1 (|2000f) to obtain the following stellar properties: M* 

0.54 
0.37 



1.36 ± 0.04 M , Age 



L* = 3.91 and log gi so = 4.19 ± 0.03. The model isochrones are shown in Figure[5] 



The value of log g obtained from the isochrone fit is somewhat larger than that found from spectral 
analysis, which was also noted by T10. We investigated the effect this higher logg would have on the fit 
to the spectral line shapes by fixing log g spec to the isochrone value and reassessing the spectral analysis. 
We found that in order to retain a good fit to the Mg I b lines, the magnesium abundance [Mg/H] must 
be reduced by 0.12 (in a similar fashion to WASP-1, see Stempels et al.l 120071 ) due to the inherent anti- 
correlation between [Mg/H] and log g. The degeneracies between the spectroscopic parameters do not rule 
out the higher log g value as long as [Mg/H] is allowed to be under-abundant, and we take the uncertainty 
in log g to include both values. Altering log g and [Mg/H] in this fashion had no significant effect on the 
derived effective temperature or metal abundance. 



3.2. Planet parameters 



To determine the properties of the planet HAT-P-14b, we simultaneously modell ed the light curves and 
radia l velocities w i th a g lobal MCMC fit. Details of this process are described in iCameron et all (|2007r ) 
and IPollacco et al.l (120081) . The free parameters in the fit are: orbital period P; transit epoch To; transit 
duration T^ UI ; squared ratio of planet radius to star radius Rp/R^; impact parameter b; RV semi-amplitude 
K; Lagrangian elements ecoso; and esinw, where uj is the longitude of periastron; and the systematic 
offset velocity 7. In this particular case, two systematic velocities were fit to allow for instrumental offsets 
between the SOPHIE and FIES datasets. These were found to be: 7sophie = -20.3834 ± 0.0027 km s _1 and 
7fies = -20.4290 ± 0.0017 kms" 1 . 

The chain length was 10,000 points and the resulting best-fit parameters and uncertainties are shown 
in Table [3] The photometric error bars were re-scaled to take into account any underestimation in the 
uncertainties so that Xred = X 2 /dof = 1 (dof = number of points - number of fitted parameters). We used 



Fig. 4. — Adaptive optics image of HAT-P-14 in the i^-band. The corrected PSF has a FWHM of 3 pixels 
(0. 12"). No contaminating sources are detected to a limiting magnitude of H ~ 17.1 mag at a distance of > 
0.48" from the centre of the corrected profile. 




proxy R*jMj against effective temperature for [M/H] = 0.08. The solid lines indicate (from left to right) 
the ages 0.01, 1, 1.6, 2.5 and 5 Gyr and the dotted lines show masses of 1.5, 1.4 and 1.3 M Q . The position 
of HAT-P-14 is indicated within the error ranges. 
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the IClarea (120001 . I2004T) non-linear limb darkening coefficients for the appropriate stellar temperature and 
photometric passband appropriate for each light curve. 



In order to estimate the stellar mass and properties derived from it , we used the Torres et al.1 (feOlOal) 
empirical calibrations from eclipsing binaries. This method, adapted by Enoch et all (|2010h . computes the 
stellar mass as a function of T e g, log and [M/H] and has the advantage of being independent of stellar 
models. We find M*= 1.30 ± 0.03 M Q , similar to the stellar mass derived from the isochrone analysis 
undertaken here (M*= 1.36 ± 0.04 M , Padova models) and in T10 (M*= 1.39 ± 0.05 M , YY models). 

We performed the MCMC fit and derived the following planetary parameters: M p = 2.25 ± 0.06 Mj, 
R p = 1.19 ± 0.05 Rj and p v = 1.33 ± 0.16 pj. In addition, we analysed the data solely from T10 and found 
that the fitted parameters were completely consistent with their results. Finally, we analysed all the datasets 
simultaneously in order to optimise the best-fit solution using all the information available. The results are 
shown in Tableland are in good agreement. 



3.2.1. Eccentricity 

In an initial fit to the SOPHIE and FIES RVs, the eccentricity was allowed to float, yielding e = 0.054 
± 0.017. We analysed the datasets independently and found that a jitter term of 15 m s -1 was required to be 
added in quadrature to the FIES RV uncertainties in order to obtain Xred = ^ • ^ n ^ s ma y ^ ue to astrophysical 
noise, as also noted by T10, and/or instrumental systematic noise which has not been accounted for. The 
FIES RVs do not significantly detect an eccentric orbit, yielding e = 0.04 ± 0.03, but equally do not rule one 
out. 

In contrast, the SOPHIE RVs do not require any jitter term to reconcile the rms scatter with the mea- 
sured uncertainties. Fitting the SOPHIE RVs alone, a more significant detection of eccentricity was found, 
e = 0.081 ± 0.022 with x 2 = 2.0 for 9 RV points and 4 fitted parameters (K, ecosu, esinw, 7sophie)- 
Fitting a purely circular orbit to the same data, we found Xchc = 22.6. Applying the Lucy-Sweeney test 



(ILucy & Sweeneylll971h . the eccentricity is detected at the 3a level. 

A joint fit to both datasets reveals e = 0.06 ll^o, X 2 = 34.4 and xl iTC = 49 for 33 RV points and 5 
fitted parameters (including two systematic velocities, 7sophie and 7fies)> indicating the eccentricity is 
non-zero at the 2.7<t level. We also apply the test to the T10 Keck RVs, including their 7.3 m s _1 jitter term. 
We find the eccentricity is detected at the 4a level, e = 0. 104 ± 0.012 , with x 2 = 12.3, x 2 iic = 83.9, 14 RVs 
and 4 fitted parameters. By combining all three datasets, our analysis suggests that the orbit is not circular 
with a 6a confidence, e = 0.095 ± 0.011 (x 2 = 55, Xcirc = 133, 47 RVs and 6 fitted parameters (including 
three systematic velocities, 7sophie, 7fies and 7 Ke ck)- 



Table 3: System parameters for HAT-P-14 



Parameter (Unit) Value Torres et al. (20 10b) Combined 
Photometric parameters: 

P(d) 4.627653 ± 0.000008 4.627669 ± 0.000005 4.627657 ± 0.000005 

T (HJD) 2455301.03431 ± 0.00036 2455134.43854 ± 0.00026 2455134.43856 ± 0.00027 

T dur (d) 0.0938 ± 0.0013 0.0912 ± 0.0017 0.0932 ± 0.0010 

AF = i?2/R2 0.00678 ± 0.00016 0.00648 ± 0.00024 0.00671 ± 0.00014 

Rp/R* 0.0824 ± 0.0010 0.0805 ± 0.0015 0.0819 ± 0.0009 

b 0.915 ± 0.005 0.891±g;{$| °- 907 ± °- 004 

a/R* 8.60 ±0.31 8.87 ±0.29 8.49 ± 0.25 

83.5 ±0.4 83.5 ±0.3 83.2 ± 0.2 
Spectroscopic parameters: l 

Kikms- 1 ) 0.2291 ±0.0051 0.2190 ± 0.0033 0.2224± 0.00261 Z 

ecosw -0.025 ±0.01 -0.009 ± 0.009 -0.0 14± 0.006 1 

esmu 0.056 ±0.025 0.106 ±0.013 0.094 ± 0.012 

e 0.061 ± 0.022 0.107 ± 0.013 0.095 ± 0.011 

w(°) U4tn 94 ±4 99 ± 4 
Derived parameters: 

M P (M Jup ) 2.25 ± 0.06 2.23 ± 0.06 2.20 ± 0.04 

Rp(R Jup ) 1.19 ±0.05 1.15 ±0.05 1.20 ± 0.04 

pp(pj) 1.33 ±0.16 1.47 ±0.20 1.28 ±0.12 

log5p(cgs) 3.56 ±0.03 3.62 ± 0.04 3.54 ± 0.03 

o(AU) 0.0591 ±0.0004 0.0606 ± 0.0007 0.0594 ± 0.0004 

r e q,A=o ( K ) 1585 ± 35 1570 ± 34 1597 ± 29 
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4. Conclusions 



We report the independent discovery and follow-up observations of the transiting planet HAT-P-14b 
found by T10. We have analysed the follow-up data used to secure the detection, and derived an independent 
set of planetary parameters which are consistent with those found by T10. We provide additional evidence, 
through AO observations and additional bisector span measurements, that the planetary signal is not caused 
by astronomical false positives. 

We find the properties of the parent star to be in good agreement with those found by T10, which is 
reassuring given the different techniques used in their derivation. The star is a hot F5 dwarf with T e g = 6583 
± 100 K and appears to be relatively you ng at 1.61° \ Gyr, w ith a slightly super-solar metallicity [M/H] = 
0.08. Using the empirical calibrations of iTorres et al.l (I2010ah we determine the stellar mass to be 1.30 ± 
0.03 M©, which agrees well with the values we derive from theoretical stellar models. 

By combining data from both groups, we are able to refine the orbital ephemeris P = 4.627652 zfc 
0.000004 d, T = 2455134.43823 ± 0.00022 d, and planetary parameters. HAT-P-14b has R p = 1.20 ± 0.05 
R.t, M p = 2 . 20 zfc 0.04 Mj and p p = 1.28 zfc 0.12 pj. It appears to blur the bi-modal distribution noted by 
Bakos et all (120 101) in which massive (M p ~ 2 Mj) planets appear either inflated 1.3 Rj) or compact 

1 Rj). We find the orbit to be significantly eccentric, e = 0.095 ± 0.01 1, which hints at the possibility 
of another body in the system. 

As noted by T10, the system promises excellent opportunities for further follow-up and characterisa- 
tion, in particular due to the brightness of the host star, V m3ig = 9.98. The star is relatively fast rotating 
for a planet host, vsini = 8.4 zfc 1.0 km s" 1 , making it a good candidate for the detection of the Rossi ter- 
McLaughlin effect (RM. lRossiterlll924j : lMcLaughlinlll924T) . However, due to the very high impact param- 
eter, b = 0.905 =fc 0.005, the RM am plit ude is only expecte d to be ~ 20 m/s. Despite this, it is still an 
attractive target as [Winn et all (120101) and lSchlaufmanl (120101) have shown that misaligned planets appear to 



preferentially orbit stars with high T c g and M± such as HAT-P-14. 

The high temperature of the host star places HAT-P-14b in the highly-irradiated 'pM' category of 
Fortney et all (|2008l) and is a good candidate for investigating temperature inversions and day/night con- 
trasts. An observation of the secondary transit would be particularly interesting to further constrain the low 
eccentricity. However, the orbital geometry of the system means that HAT-P-14b is on the borderline of 
having no secondary eclipse. The results from the joint analysis predict that it will not occur with more than 
3cr confidence. However, the fit to the data presented in this paper alone suggests that the secondary eclipse 
could be grazing. Since these opposing secondary eclipse predictions arise from extremely small differences 
in orbital eccentricity, whether or not we observe such an transit will provide further constraints on the sys- 
tem geometry. HAT-P-14 thus provides excellent opportunities for future investigations and informative 
characterisation. 
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